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SUMMARY 


An examination of oscillatory stability for a straight-winged airplane 
with large concentrated wing-tip masses was made using wing-bending and 
airplane-pitching degrees of freedom and considering only quasi -steady 
aerodynamic forces. It was found that instability caused by coupling of 
airplane pitching and wing bending occurred for large ratios of effective 
wing- tip mass to total airplane mass and for coupled wing-bending frequen- 
cies near or below the uncoupled pitching frequency. Boundaries for this 
instability are given in terms of two quantities: (l) the ratio of effec- 

tive tip mass to airplane mass, which can be estimated, and (2) the ratio 
of the coupled bending frequency to the uncoupled pitch frequency, which 
can be measured in flight. These boundaries are presented for various 
values of several airplane parameters. 


INTRODUCTION 


The placement of fuel or armament in wing-tip pods is a practice which 
has been in existence for some time on fighter-type aircraft, and there are 
several ways In which large concentrated masses might appear at the wing 
tips of future aircraft. There is the possibility that in nuclear powered 
aircraft, where the personnel and power plant must be well separated, this 
technique will be employed. Long range aircraft with supersonic dash capa- 
bility might carry disposable wing-tip sections with large fuel tanks. 

Also the mounting of engines at the wing tips is proposed for certain air- 
plane types such as vertical take-off airplanes. 

The natural wing-bending frequency of an airplane might be lowered by 
the presence of large wing-tip masses to the extent that it is near or 
even below the natural pitching frequency of the airplane under certain 
flight conditions. Under these conditions the effects of coupling between 
pitching and wing bending are of concern, since the damping ratios of 
structural modes are inherently low and any unfavorable coupling might 
cause oscillatory instability in coupled wing bending. It is the purpose 
of this study to establish conditions under which such wing-bending 
instability is likely to occur. 
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The analysis is made by means of equations of motion with airplane 
pitch and fundamental wing-bending degrees of freedom. Unsteady aerody- 
namics, that is, the effect of frequency on aerodynamic coefficients, and 
structural damping are neglected. 
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NOTATION 


derivative with respect to real time, — 

dt 


airplane pitching moment of inertia about the center of gravity , 


slug ft 
total airplane mass, slugs 


rate of change of pitching moment with pitch angle, 0, 
f' 1 1 y f t — lb/radi an 


rate of change of pitching moment with pitch velocity, 0, 
ft -lb sec/radian 

generalized aerodynamic pitching moment, ft -lb 
generalized aerodynamic wing-bending force, lb 
wing area, sq ft 

forward velocity of airplane, ft/sec 

rate of change with the angle ^ of the generalized wing-bending 
force due to an angle -of -attack disl ribution ^ a Q (r|), 


2 / l a 0 a o^- T 1^ lb/radian 


*a Q 

Z0 


rate of change of generalized wing-bending force with pitch 


angle 0, 2 P 2ga 0 dT], lb/radian 
’-'o 


Ye 

Z0 
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rate of change with the angle — of the generalized vertical 

' * 

force due to an angle-of-attack distribution — a Q ( rj ) , 

2 / Z a dTj, lb/radian 
'-’o u 


rate of change of generalized vertical force or lift with pitch 

angle 9,2 Z^dr^ lb/radian 
^o 

spanwise bending deflection as a function of t] referenced to 
node and normalized on cantilever wing-tip deflection y 

spanwise bending deflection as a function of tj referenced to 
wing root and normalized on cantilever wing-tip deflection y, 
a(n) - a(0) 

value of a(r)) at T] = 0 or the wing root 
mean aerodynamic chord 
Naperian base 


reduced pitching frequency parameter, — rj— 


rate of change with the angle — of the spanwise aerodynamic 
load distribution due to an angle-of-attack distribution 
“ a ( t] ) , lb /radian 

y 

rate of change with the angle ^ of the spanwise aerodynamic 
load distribution due to an angle-of-attack distribution 


y a 0 (n), lb/radian 


rate of change with Q of the spanwise aerodynamic load distribu- 
tion due to a pitch angle 0, lb/radian 
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effective wing-tip mass (both tips), 2 / pa 0 dTi, slugs 


generalized wing-bending mass, 2 / pa 2 dr), slugs 


derivative with respect to dimensionless time. 


dynamic pressure, lb/sq ft 


pitching radius of gyration 


/I 

n/M 9 


spanwise distribution of static mass moment about airplane center 
of gravity, slug ft 

real time, sec 


static margin, — , ft 
In 


longitudinal distance from airplane center of gravity, positive 
forward, ft 

wing-fuselage local aerodynamic center relative to airplane 
center of gravity, positive forwarc. 

wing local aerodynamic center relative to airplane center of 
gravity, positive forward 


longitudinal tip mass center of gravity relative to airplane 
center of gravity, positive forwarc , ft 
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y 

n 

9 



n(n) 

T 


a 


wing-tip deflection relative to wing root (cantilever ) , positive 
down, ft 

damping ratio of coupled wing-bending mode 

damping ratio of uncoupled, free-free, wing-bending mode 

damping ratio of uncoupled pitching mode 

dimensionless spanwise station 

pitch coordinate, positive nose up, radian 

component of pitch in phase with bending velocity during undamped 
oscillation 

spanwise mass distribution of airplane, slugs 
dimensionless time, cjgt 

WB 


u>B 

U)y 


undamped natural frequency of coupled wing-bending mode, 
radians/ sec 

undamped natural frequency of uncoupled, free-free, wing-bending 
mode , radians / sec 

undamped natural frequency of uncoupled pitching mode, J — j“ , 
radians/sec 


Dots over symbols are used to indicate differentiation with respect 
to real time. 


DERIVATION OF EQUATIONS OF MOTION 
Selection of Degrees of Freedom 


A system with two degrees of freedom is considered for the purpose 
of studying oscillatory instability caused by coupling of wing bending 
with airplane-pitching motion. The selection of the two degrees of free- 
dom was based on the results of some four degree of freedom studies on 
specific example airplanes, where the degrees of freedom were airplane 
pitch, airplane vertical translation, cantilever wing bending, and canti- 
lever wing torsion. The examples were similar in character to the type 
of airplane considered in the present analysis. The development of the 
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equations was by the method shown in appendix C of reference 1. In these 
studies conditions of undamped oscillation in inodes which were primarily 
wing bending were found and the modes of oscillation were examined. 

Two important facts were noticed: 

1. The relation between vertical translation and cantilever 
wing bending was nearly the same as would be expected for the 
airplane in an uncoupled, free-free, wing fundamental bending 
mode ( in vacuo ) . 

2. The role of wing torsion in reducing the generalized 
coupled mode damping force to zero was found to be small in 
comparison with the role of pitch, and instability could be 
achieved without the inclusion of the torsion degree of freedom 
under conditions not greatly different than with it. 

The above observations indicated that the essential mechanism involved 
in this oscillatory instability could be preserved in a system of only two 
degrees of freedom, airplane pitch and a wing-bending degree of freedom 
based on the uncoupled, free-free, wing f irsn-bending mode in the absence 
of aerodynamic forces. Such a simplification is of great help in the 
examination of the general nature of the instability. 


The Generalized Wing-Bending Coordinate 


Before the generalized bending coordinate can be established, a brief 
look at the nature of the fundamental, uncoupled, free-free, wing-bending 
mode (in vacuo) is necessary. Consider an airplane with large concentrated 
wing-tip masses oscillating in such a mode. If the generalized wing- 
bending coordinate at any time is called y, and the corresponding mode 
shape is called a(r)), the vertical position of the wing relative to the 
node at any span station t] is ya(rj) (fig. 1). The mode shape a(r]) may 
be considered to be made up of two parts 

a (n) = a-oCl) + a (o) 

where a 0 (rj ) is simply the mode shape referenced to the wing root rather 
than the node (fig. 2). With this definition and with a 0 (r|) normalized 
to 1.0 at the wing tip, it is evident that j is the cantilever wing-tip 
deflection (see fig. l). The quantity a(0) can be determined from the 
equilibrium requirement for the freely oscillating system in the absence 
of aerodynamic forces 





n(T]) a (Tj)dT] 


0 
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where p(t] ) is the spanwise mass distribution including the concentrated 
masses. Because the mode is symmetric the integration need only be over 
the semispan. Substituting a(ri) = a o (r)) + a(0) and dividing the equation 
by y yield 


or 


ap dr] = 2 / a G p dt] + 


2a(0) J p 


dn 


0 


a G p d^ + a(o)M = 0 


where M is the mass of the airplane. Then 


where 


a(ri) = a 0 (ri ) 


m 

M 


m 


2 



l 

a c p dt} 


(1) 


It should be noticed that if the distributed mass in the wings were con- 
sidered negligible, m would simply be the sum of both concentrated wing- 
tip masses. (Note in fig. 2 that ag(o) = 0, excluding the fuselage 
mass, and a (l) = 1.) Now if a 0 (Ti) is assumed to remain unchanged, 
a(n) changes with m/M only as shown in equation (l), greatly simplify- 
ing the use of a free-free mode as a degree of freedom. Fortunately the 
present analysis is fairly insensitive to small changes in a 0 (q), and 
this assumption has little effect on the results. The reason for this 
will become evident later. 

The generalized wing-bending coordinate, then, will be represented 
by the cantilever wing-tip deflection y, where a value of y implies a 
spanwise deflection relative to the node of y (n 0 (T)) - 


Kinetic Energy 


(i/ 2 >/ 


The kinetic energy of the two degree of freedom system 


is 


v 2 dM where v is the local velocity in space at the location 


of the” differential mass dM and the integration is made over the whole 
airplane. Small angles being assumed, the local velocity may be expressed 
in terms of the two coordinates and can be expressed as 


v = 0x - ya(-n) 
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where x is the longitudinal distance from the airplane center of gravity, 
positive forward. Then 


KE = | 


0 2 x 2 - 20xya(i]) +y 2 a(T|) 2 &M 


airplane 


d 


= ^ 0 2 I - 2y0 P sa di) + y 2 P |j.a 2 d.r] 


1 

2 


where I is total pitching moment of inert:. a about the airplane center 
of gravity and s(r)) is the spanwise distribution of static mass moment 

about the airplane center of gravity. Note that / sa drj = / sa 0 dT] 

p i Jq Jo 

since / s dr\ = 0. Thus, letting 


1 


Ga Q dr) 


x p - 


f 


pa 0 dr) 


mg = 2 / a 2 jj. d*Q 

Jq 


and remembering that m = 2 / a. Q \x dr) 


KE = 7j 0 2 I - y0Xpm + \ y 2 ^ 


When the distributed mass of the wing is considered negligible, Xp 
simply becomes the longitudinal distance between the center of gravity 
of the concentrated wing-tip mass and the airplane center of gravity. 


Potential Energy 


The only potential energy in the systen is the strain energy in the 
wing. Since ao(r)) is held constant the strain energy is a function of 
y only and can be written 

< 

PE - | Kyy 2 
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where K y is an effective spring constant. In terms of the uncoupled, 
free-free, natural bending frequency w y , K y can be written as rrigtOy . 
The bending shape is very critical in the determination of Kyj but, as 
will be seen, this problem is bypassed in the method to be presented. 


Lagrange ' s Equations 


The two simultaneous equations of motion are obtained from Lagrange's 
equation by substituting the expressions for kinetic and potential energy 
with the coordinates y and 0. 


_d_ Ake'N Bra 

at \Be/ Be 


d /Bke\ Bra „ 

dt + By = 


I 


(2) 


where BPE/Be happens to be zero. The Q's are the generalized aero- 
dynamic forces. Structural damping is neglected. When the indicated 
differentiations are performed, the equations become 

10 - mxpy = ( 3 a ) 

-mxp0 + m g y + K y y = Q 2 ( 3*> ) 


Generalized Aerodynamic Forces 


A generalized force is the work done per unit displacement when the 
system undergoes a virtual displacement of one of the degrees of freedom. 
The generalized aerodynamic force is assumed to be composed of terms 
which are linearly dependent on y, 0, and their derivatives with respect 
to time. If unsteady aerodynamics are neglected, the essential components 
of the generalized forces for a straight -winged airplane can be written 


1 l 


a 


Qi = Mg0 +M e 0 - l 2 / — x ac dri hr 

J 0 


Qp = 2 


al e dr^) 


0 + 12/ a — dT) Jy 


where a term is positive if the force or moment involved is directed in 
the positive direction of the appropriate coordinate. The derivatives 
M 0 and M 0 can be readily written in terms of the conventional airplane 
stability derivatives . 
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M 9 = iScCCmq + C^) , % = qScC^ 

The term l a (r\) can be interpreted as the air load distribution resulting 
from a spanwise angle-of -attack distribution defined by a(rj) in radians, 
and Iq (t) ) can be interpreted as the air load distribution resulting from 
a uniform angle-of-attack distribution of 1 radian. The generalized bend- 
ing force due to 9 and pitching moment due to y (from downwash lag) 
were considered negligible. 

The air load distribution l a can be written as the combination of 
the loading caused by an angle-of-attack distribution of a 0 (T]) radians 
and (m/M)lg or 


= ia D - 1 19 


m 


The advantage of this approach is that the generalized aerodynamic forces 
can be written as explicit functions of m/M. If this relation and the 
fact that a(i]) = a 0 (Tj) - (m/M) are used, the expressions of Q x and 
may be rewritten 


Qi = Mq0 +M00 - i (x a 2 J l ao dTi - g 2 J l e >. ac &n j 


Q 2 = 2 


m 


Z 0 dT l " M 2 / l 9 dT \ 9 + 


1 

V 


2 



a o 2 a 0 dl l " M 




i 

l e dr) 


y 


where the constant x a in Q 1 results from the fact that Z aQ contains 
no tail or fuselage effect and the assumption that x a is uniform over 
the span. For convenience all the terms involving spanwise integrations 
will be designated as aerodynamic derivatives represented by single 
symbols with subscripts. 


Qjl = 


1 

V 


x a z a r 


IT1 w 

+ M % 


(4a) 


q 2 



a 7 

m z e 


9 + 


i 

V 




(4b) 


The term Z q is the rate of change of vertical force with 9 and 
can be written in terms of the familiar airplane lift curve slope as 
Zg = -qSCL a ; Za Q is the rate of change of vertical force caused by an 
angle-of-attack distribution (y/v)a 0 (Tj) with y/V; Y g and Y ao 


are 
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generalized wing-bending force derivatives whose definitions follow from 
the subscripts. Small changes in a 0 (r|) have been found to have little 
effect on these derivatives. 


Equations of Motion 


If equations ( 3 ) and (4) are combined a set of homogeneous equations 
of motion can be written. With D as the differential operator d/dt, 
these equations are 


( ID 2 


MgD - M Q JG + 


I-. 2 1 ( m 

_-mXpD 2 


V Vm 


a^a. 


'D 


y = 0 

(5a) 


-mXpD -lY Q 


m 

M 


0+ im 


igD 2 - 


1 

V 


m 


a G 


- Ml y e +Z 


a G 


§> z < 


D + Ky ("y = 0 


(5b) 


These equations will be useful in the examination of the wing-bending 
mode but, since no independent vertical translation of the airplane is 
allowed, the airplane short-period mode is not well represented. In the 
development of the equations the coordinate deflections are held to small 
displacements, unsteady aerodynamics are not included, and a constant 
bending shape a 0 (r|) is assumed. 


DETERMINATION OF STABILITY BOUNDARIES 
Dimensionless Form of Equations of Motion 


The reduction of the equations of motion to the two degree of freedom 
form of equations ( 5 ) makes possible an analysis in general algebraic 
terms of any airplane to which the equations are applicable. For a 
general study it is convenient to group the variable quantities describ- 
ing the airplane characteristics and flight conditions into dimensionless 
parameters . 

In the analysis of equations ( 5 ) the following dimensionless quantities 
were found to be convenient: 

r dimensionless time, Ugt 

m’ effective tip mass ratio, ^ 

ru)Q 

k 0 dimensionless frequency in uncoupled pitch, 



12 


1 1 M£ 

damping ratio of airplane in uncoupled jitch, — 

xl dimensionless effective tip mass center- of -gravity location, ^ 
Jr r 

. x a 

x^ dimensionless wing center-of-pressure location, — 

Me 


u* dimensionless stability margin, 

mg 

mg generalized wing mass ratio, — 


rZe 


ZI 

a o Z, 


a 0 


e 


Y a z e 


Ye 

Ya r 


Y 1 

a o z e 


In terms of these new quantities equatiors (5) become 


P 2 e + 2£ 0 d0 + e - k^m- + ^ z a 0 )p(?) = 0 

(6a) 


nrx-J 


m 


+ ^rrr («- ■>•)» + 


g 


u ' m 


ke 

+ miu 1 


m 1 


•l^+Z- ) + m 


,t2 


iV^Yf, = 0 

(6b) 


J P W + W' 


where p is a new differential operator meaning d/dT. 


The Criterion for Neutral Stability 


The conventional method of determining the boundaries of neutral 
stability is by the Routhe criterion. However, even with a fairly simple 
fourth-order system such as this the computations are quite tedious and 
any physical feeling for the mechanism of the instability is lost. The 
equation system (6) lends itself very well to another method of determin- 
ing stability boundaries which is satisfactory for the present analysis. 
It greatly aids in physical interpretation of the cause of instability, 
and the computations required are very simple. 
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The system is assumed to he in a steady -undamped oscillation at the 
coupled natural wing-bending frequency, and the necessary conditions for 
this undamped oscillation to exist are used to define the stability 
boundaries. If the coupled natural frequency in wing bending is Wg, 
then based on the dimensionless time used in equation (6), the natural 
bending frequency is ojg /^q> which will be called ft. 

By assuming an undamped solution to equations (6) of the form 
y/r = (y 0 /r)e iaT and substituting this solution into the pitching-moment 
equation (6a) it is easily shown that 


e = (A+iB) ^ e iftT 


where 


A = 


m'x£ft 2 (l - ft 2 ) + 2k 0 (^rn' + -jf 1 


(1 - ft 2 ) 2 + (2 £ 0 ft)‘ 


and. 


k e( m ' + TT Z a Q )^ 1 " ft 2 ) +2m'x^ 0 ft 3 
B_ (1 - ft 2 ) 2 + ( 2 £ 0 ft ) 2 

But since y/r = (y 0 /r)e iaT and p(y/r) = ift(y 0 /r )e iaT , Q can be written 


and 



P 2 e = -ft 2 e 

Using the expression for 9 , the generalized bending force equation (6b) 
can de written as 



B(Y^ - m» ) 


+ 


k e 

m gU' 



X- 


mi 


Aft 2 


A 

u'mL 


(Y> - m* ) + 


Ay 


w j 


0 


(7) 
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Since this equation is valid only for the case of undamped oscilla- 
tions (neutral dynamic stability), the coefficient of p(y/r) must vanish. 
Thus 




B(Y£ -m») 
ftu* 


+ 


^9 

u’ 



( 8 a) 


This leaves, in equation (7), for sinusoidal motion 


x£m 


tut 


m; 


u'rni 


(Y* - m* ) + ( 7 — 




= 


( 8 b) 


Equation ( 8 a) defines the stability boundary in terms of 0. rather than 
Wy/w 0 . For design use, Wy/wg could easily be found in terms of Q from 
equation ( 8 b) provided is known accurately. Only equation ( 8 a) was 

used in this study and it should be noted that neither m£ nor Wy appears 
in this equation. Hence, Ky does not enter in any fashion and the stiff- 
ness problem is bypassed, as mentioned earlier. Such an approach, which 
does not employ equation ( 8 b) to determine Uy/w 0 , would be useful if the 
analysis were used in conjunction with actual flight test, where the 
coupled bending frequency ojg (wg = au Q ) would be measured. 


Stability Boundaries 


Equation ( 8 a) has been used to plot stability boundaries on a plane 
of m’ versus ft with all other parameters held constant (figs. 3 , 4 , 
and 5). Since the number of independent parameters is large, presenta- 
tion of boundaries for all possible combinations of these parameters over 
a range of values would be prohibitive. For this reason only a few combi- 
nations are shown to illustrate the general effects of variations in each 
quantity on the boundary. 

In figure 3 the quantities varied are x* k 0 , and-u'. For computa- 
tion of these boundaries the other parameters were held constant at the 
following values: = 0-35, = 0, Y£ = 0.270, = O. 255 , = 0.108. 

All these except x^ are values which were calculated for an example air- 
plane and thought to be typical of a straight- rfing airplane with large 
wing-tip pods. The value of x^ was chosen partly because of simplifica- 
tions which accompany the value, and partly because it is well within the 
range of reasonable values. Effects of variations in these quantities 
are dealt with later. The values of k 0 and a* were chosen to fall in 
what was considered the range of interest for airplanes which are most 
likely to encounter this problem. 

The boundaries in figures 4 and 5 are calculated with x£ = 0, 
mainly because of the great simplification afforded by this value. When 
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is zero, the boundaries become independent of k^ and u ! as is seen 
later. The significance of the range of values of x a /u used in fig- 
ure 5 will be discussed in the next section. 

DISCUSSION OF STABILITY BOUND ABIES 
Mode of Instability 

The stability boundaries indicate conditions under which undamped 
oscillations of the system will occur. Since the system is regarded as 
having two degrees of freedom, it must be determined whether the neutrally 
stable mode is what has been called coupled wing bending or a coupled mode 
which is predominantly pitching. The nature of the neutrally stable mode 
can be determined quite easily for any particular case by substituting the 
appropriate natural frequency back into the equations and determining the 
relation between the pitch and wing-bending coordinates. 

It was found that the modes of oscillation associated with the 
boundaries considered in this study involved a maximum value of the ratio 
r|0|/|y| which was on the order of kg or, in this study, Chi . Physi- 
cally, this means that a point on the airplane fuselage a longitudinal 
distance from the airplane center of gravity equal to the radius of 
gyration would have linear oscillation amplitude with respect to the 
center of gravity, caused by pitching, usually less than one-tenth the 
oscillation amplitude of the wing tip referred to the center of gravity. 
From this observation it can be said that the coupled mode associated 
with the boundaries is primarily one of wing bending. (The largest pitch- 
ing amplitudes occur under conditions of large longitudinal displacements 
of the tip mass center of gravity from the airplane center of gravity and 
the higher frequency ratios.) 

The phasing of pitch with respect to cantilever bending deflection 
in these neutrally stable modes (except for one boundary discussed later) 
is such that 6 leads y by an amount ranging from a fairly small angle 
to something greater than 90° (vector plot, sketch (a)). The larger lead 
angles are associated with small values 
of ft and the smaller lead angles with 
values of ft near 1.0. Physically this 
phase relation means that when the wing 
tip is traveling downward with maximum 
velocity during the undamped oscillation 
(at zero displacement), the pitch angle 
Is positive. In other words, there is a 
positive component of pitch In phase with 
bending velocity. 

It is this in-phase component of 
pitch which is the major factor in 
reducing the damping to zero, since the 

Sketch (a) 
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air loads caused by the in-phase pitch effectively alter those which 
normally supply the damping in the uncoupled wing-bending mode. This 
fact can be illustrated with the aid of a rather simplified picture of 
the air loads which provide the damping. The distribution of vertical 
velocity along the span, ya(r}), for an airplane with rather large wing- 
tip masses might appear as in figure 6(a) during a wing-bending oscilla- 
tion. Such a velocity distribution would have associated with it a 
spanwise angle -of -attack distribution equal to (y/V)a(r]) (fig. 6(b)) during 
forward flight and consequently a spanwise loading (y/v)Z a which would 
resemble that in figure 6(c). In an uncoupled bending oscillation or a 
bending oscillation in which there were no pitch in phase with bending 
velocity (Oy) , it is this loading which would tend to damp the oscillation. 
It can be seen from the figures that such an oscillation would always be 
stable since the air loads oppose the motion essentially over the entire 
span. 


In the unstable bending modes mentioned above, however, there is a 
component of airplane pitch in phase with bending velocity which is 
accompanied, of course, by a spanwise loading similar to that shown in 
figure 6(d). The effective span loading in phase with y, which deter- 
mines the damping, then, becomes the sum of the loadings of figures 6(c) 
and 6(d) shown in figure 6(e). It is easily seen from a comparison of 
figure 6(e) with 6(a) that the air load favors rather than opposes the 
direction of motion of the wing over the center part of the span. Integra- 
tion of the product of the loading of figure 6(e) with the mode shape 
a(r|) over the span indicates whether the neu effect is stable, neutrally 
stable, or unstable oscillation. Again by looking at the figures one can 
see that the stability of the oscillation as. determined from this integra- 
tion will depend mostly on the amount of pinch in phase with y and the 
spanwise location of the nodes. This latter quantity is determined by the 
effective tip mass ratio m* and the bending shape. Movement of the nodes 
toward the tips (or increasing m* ) and increasing the component of pitch 
in phase with bending velocity are both contributory to instability in 
this mode. 

A brief look at the stability criterion (eq. (8a)) will help in the 
physical interpretation of the effects of the various parameters on the 
stability boundaries. The third term of th:.s expression is simply the 
discriminant or stability criterion for the uncoupled wing-bending mode 
(see coefficient of p(y/r) in eq. 6(b)). There is no positive value of 
m T between 0 and 1 for which the bracketed expression becomes 0 with 
realistic aerodynamic derivatives, meaning nhat the single degree of free- 
dom wing-bending system cannot become unstable. This was demonstrated 
physically above. 

The second term of equation (8a) is proportional to the component of 
pitch in phase with bending velocity since it was shown that 
9 = A(y/r) + (B/n)p(y/r) and is the only additional term which would 
appear in the stability criterion of the coupled system if x£ were 
zero. Also, it was found in the calculation of the stability boundaries 
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that the first term, proportional to pitching acceleration in phase with 
bending velocity, is not an important contributor for ft less than 1.0 
even when xi is not zero, so that the examination of only the second 
and third terms can give much insight into the behavior of the boundaries 
at the lower frequency ratios. If the first term is neglected equa- 
tion (8a) can be rewritten in the following manner: 






+ Z L 


+ 



+ m* 


2 


= 0 


(9) 


The term B/ftkg can be shown equal to 0y/ (y/V ) , where 0y — (B/ft)p(y/r) 
is the component of pitch in phase with bending velocity and y/V is 
simply the characteristic angle of attack caused by wing-bending velocity. 
As 9y goes to zero equation (9) becomes the same as the third term of 
equation (8a). The effect of this term is to alter the aerodynamic 
derivatives as they appear in the third term of equation (8a). 

A qualitative plot of m' required for neutral stability versus 
0-/(y/v) based on the simplified criterion, equation (9), is shown in 
figure 7. Since all the boundaries in figures 3, b, and 5 except one 
are associated with modes of oscillation involving pitch in phase with 
bending velocity the right-hand boundary will be of most interest. It 
can be seen that as the in-phase component of pitch is increased from 
zero, a value is reached where neutral stability can exist for O^m' <C1.0. 
As 0y/(y/v) increases beyond this value, the value of m' required for 
neutral stability decreases. This figure will be referred to as the 
effects of the various parameters are discussed. 

The general effect of the first term of equation (8a) is to raise 
the boundaries in figure 3 with positive values of x£ above what would 
be found with equation (9) in the region of ft = 1.0 and above. 


Relation of Boundaries to Existing Airplanes 


Now that stability boundaries have been established in terms of 
several dimensionless quantities, it is of interest to know where exist- 
ing airplane types stand relative to these boundaries. First, orienta- 
tion on the m' - ft plane will be considered. For airplanes with no mass 
concentrations at the tips, the effective wing-tip mass is quite small 
compared with the total airplane mass (m 1 probably less than 0.1 ), and 
the fundamental wing-bending frequency is high in comparison with the 
uncoupled pitching frequency of the airplane. Such airplanes would be 
below and to the right of the portion of the m' - ft plane shown in fig- 
ures 3 and One possible exception to this might be a flying wing 
configuration where the total mass is distributed along the wing. In 
this case m' might be fairly large. (The value of m’ for a uniform 
beam, an extreme example, is about O. 38 .) 
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The addition of concentrated masses to the wing tips tends to move 
a given airplane upward and to the left on the plot of m’ versus SI, and, 
as the concentrated wing-tip mass is increased, the trend continues in 
the same direction toward the unstable regions. However, even airplanes 
which are now considered to have very large wing-tip pods have insuffi- 
cient values of effective tip mass ratio m' to reach any of the bounda- 
ries shown in figures 3 and 4. To cite a specific example, a fully loaded 
F- 89 D, which represents nearly the upper limit in wing-tip mass ratios for 
current airplanes, could be shown as a point slightly above m* equal to 
0.25 and near Si = 1.0 for a condition of high speed and low altitude. 

The location of the boundary itself on the m' - SI plane is a func- 
tion of some quantities which should be related to existing airplanes. 

The quantities x£ and u' are fairly straightforward. They are dimen- 
sionless on the radius of gyration instead of the more familiar mean aero- 
dynamic chord, but these two characteristic lengths are of similar size. 
Again using the F- 89 D as an example, the approximate range of Xp is 
from 0.25 to -0.4, and of u' is from nearly 0 to 0.45- The large value 
of u* is the result of a center-of -pressure shift in the transonic 
range, and similar values can be expected for airplanes flying 
supersonically. 

The term kg is comparable in value to the familiar reduced fre- 
quency encountered in the study of classical flutter, the only difference 
being the use of radius of gyration instead of semichord as the charac- 
teristic length. A value of kg equal to 0.1 represents almost the 
upper limit for present-day fighter-type airplanes. The highest values 
of kg would be encountered at low supersonic speeds and low altitudes. 

The term x^ is the distance of the wing aerodynamic center from 
the airplane center of gravity in fraction of the radius of gyration r. 
This quantity may be either positive or negative but will usually be less 
in absolute magnitude than u' for convention;)! airplanes. 


Effects of Tip Mass Center-of -Gravity Movement 


Figure 3 indicates the effect of tip mass center-of-gravity movement 
on the location of the boundary. As the tip mass center of gravity is 
moved forward of the airplane center of gravity from x* = 0 to x' = 0 . 5 , 
the unstable region expands downward and to the right on the m' - SI 
plane. This is a destabilizing effect in the sense that an airplane on 
the stable side of the boundary could become unstable with forward move- 
ment of the tip mass. Movement of the tip mass center of gravity in a 
negative , or aft, direction from the airplane center of gravity has the 
opposite effect. However, the danger of reaching the boundary for nega- 
tive x£, with an actual airplane would be very small, so it is not con- 
sidered. The marked effectiveness of Xp in moving the boundary suggests 
that movement of the tip mass center of gravity would be a useful tool for 
eliminating any problem involving this type of instability. 
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The effect of forward displacement of the tip mass on the mode of 
oscillation is to increase the amplitude and lead angle of 9 with 
respect to y by imposing a pitching moment in phase with bending 
acceleration. This also means that over most of the frequency range 
forward movement of the tip mass increases the amount of pitch in phase 
with bending velocity. By referring to figure 7 it can be seen that an 
increase in 9f tends to lower the value of m' required for neutral 
stability. This effect is also evident from the physical argument of 
figure 6 in that a larger loading near the midspan permits the node lines 
to be farther from the tips for neutral stability. 


Effects of kg and u' 


Unfortunately , the physical significance of independently varying 
kg and u' is not obvious. The term kg is one which arises in the non- 
dimensionalization of time and of forces and moments caused by wing-bending 
velocity. It is the key quantity in determining the local angle of attack 
at a point on the wing resulting from a given vertical velocity of that 
point. When the definition of kg is rewritten it acquires an interest- 
ing meaning which, although somewhat removed from its origin, helps in 
obtaining a physical feeling for the term. 

_ _ / Me 

Q V 4 MV 2 

The ratio under the radical in the right-hand expression may be thought 
of as a ratio of energies. Except for a factor of 2, the denominator is 
simply the kinetic energy of the airplane due to forward motion, and the 
numerator can be interpreted as the potential energy resulting from a 
unit deflection in pitch or the energy contained in an uncoupled, undamped 
pitching oscillation of unit amplitude. Thus the term becomes a measure 
of the energy contained in the uncoupled pitching mode as compared to the 
kinetic energy of the airplane in forward motion. 

The definition of u' is very simple but the way in which it enters 
the problem independent of kg is not at all obvious. Just as kg, u 1 
can be rewritten to have a physical interpretation quite apart from that 
of a dimensionless static stability margin. The term u' can be written 

M(Jg 2 r 


The numerator might be called a characteristics inertia force which would 
occur during an undamped, uncoupled pitching oscillation of unit amplitude, 
and the denominator is the corresponding characteristic aerodynamic force 
on which all aerodynamic forces have been nondimensionalized. With this 
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interpretation u' becomes a measure of the relative importance of inertia 
and aerodynamic forces. For example, an independent decrease in u 1 indi- 
cates an increase in the significance of aerodynamic forces in the system. 
(This can also be seen in eqs. (6).) 

The first point of interest is that when x£ = 0, the boundary of 
figure 3 is independent of kg and u * . The reason is that under this con- 
dition the coupling between pitch and wing bending is purely aerodynamic. 

In the more general case when x£ is not zero, both inertial and aero- 
dynamic coupling occur. 

In the range of values considered for > £ and u', kg appeared to 
have only small effect on the location of the boundaries. The small 
change that was noticed indicates that an independent increase in kg 
(u’ constant) tends to move the boundary upward and to the left on the 
m' -ft plane for positive values of x£, a favorable direction. It 
should be remembered that this observation is based only on small values 
of kg. For a given airplane kg can be changed significantly only 
through Mach number effects or change of altitude (or automatic controls). 

The effect of changes in u' on the boundaries is seen to be par- 
ticularly noticeable for values of 2 above 1.0. An independent increase 
in u’ tends to shrink the boundary to the left for positive xl. For a 
given airplane u 1 can be varied significantly only through Mach number 
effects (or automatic controls). 

Both kg and u' are quantities which are determined by considera- 
tions other than wing-bending stability, and for this reason could not 
generally be altered for the purpose of improving the wing-bending 
characteristics. 


Effects of Pitch Damping 


It should be remembered that the boundaries in figure 3 were 
calculated with £g and x a /u , and with the aerodynamic derivatives held 
constant. At this point some of the effects on the boundaries of devia- 
tions of £g and x a from their assigned values will be examined. These 
effects will be observed quantitatively only lor the condition when 
x p = 0> since under this condition the boundaries become independent of 
kg and u' , and the stability criterion becomes exactly that of equation (9). 
These results can be extended qualitatively to other conditions. 

The effect on the boundary of changes in pitch damping is illustrated 
in figure 4. The boundary for £g = 0.35 is the same as that for x£ = 0 
in figure 3, and the second boundary differs only in that £g has been 
lowered to 0.2. The effect of lowering £g Ls to extend the unstable 
region slightly as D approaches 1.0. Physically the effect of lowering 
£g is to increase the magnitude of pitch in the coupled mode near 
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ft = 1.0 through the phenomenon of resonance. This in turn increases the 
component of pitch in phase with bending velocity and hence the boundary 
is lowered. This general effect is also present where x£ is not equal 
to zero. 


Effects of Wing Aerodynamic Center Location 


The longitudinal position of the wing aerodynamic center x a enters 
into the stability criterion (eq. (8a) ) as a ratio x a /u. Aerodynamically 
the static stability term Mg is composed of two major parts, the contri- 
bution of the horizontal tail and the contribution of the wing. The wing 
contribution can be either stabilizing or destabilizing depending on 
whether the wing aerodynamic center is aft or forward of the airplane 
center of gravity. The ratio x a /u is a measure of the contribution of 
the wing to the total longitudinal static stability M 0 . A positive, 
value indicates a destabilizing wing contribution; a negative value indi 
cates a stabilizing wing contribution; and a value of -1.0 indicates that 
all the stability is contributed by the wing and none by the tail. 

The effect of changing the ratio x a /u on the boundary for x£ = 0 
is shown in figure 5. Except for the changes in x a /u, all conditions 
are the same as for the x£ = 0 boundary in figure 3. The boundaries 
for x a /u = 0.5, 0, and -0.5 give an indication of the effect of this, 
parameter for conventional airplanes. If the wing gives a destabilizing 
contribution to the longitudinal static stability, the boundaries of fig 
ure 3 are generally lowered, and the reverse is true for a stabilizing 
wing contribution. This effect is more pronounced for low values of Q 
than for high values, and, in fact, this is the only parameter discussed 
so far that has varied the location of the boundary at Q = 0 . This will 
be discussed in the next section. 

The actual mechanics of this effect are based on the fact that the 
tail contribution to static longitudinal stability is most important in 
producing a pitching moment proportional to and in phase with bending 
velocity, since the wing contribution tends to be canceled out in the 
wing-bending oscillation. The increase in 0y caused by an increase in 
pitching moment in phase with y, then, lowers the boundary in accordance 
with the effect shown in figure J. 

As a ma tter of interest a boundary was computed for x a /u = -2.0. 
This is a case where twice the total amount of static stability is con- 
tributed by the wing, or, in other words, the "tail" has a negative 
contribution as in a canard configuration. The result of this computa- 
tion was the lower boundary in figure 5. The mode of oscillation in 
this case was primarily wing bending as with the other cases, but the 
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Sketch (b) 


pitching motion lagged wing bending, 
producing a pitch component in 
antiphase with bending velocity 
(sketch (b) ) . 

Physically, the canard surface 
produces a pitching moment in anti- 
phase with bending velocity as opposed 
to the in-phase pitching moment pro** 
duced by the conventional tail. As 
a result of this, a component of pitch 
in antiphase with wing bending is 
produced which is large enough to 
reach the left-hand boundary of fig- 
ure 7* Also, by reversing the sign 
of the span loading caused by pitch 
in phase with y in figures 6(d) 


and 6(e), it is evident that zero damping can be achieved if the nodes 
are moved inward toward the fuselage (or m f reduced). The possibility 
of reaching this boundary with an actual airplane is remote because of 


the unlikely combination of low tip mass ratios 
ratios ft. 


nr 


and low frequency 


Location of Boundaries for Small ft 


For all the boundaries in figures 3 and 4 the value of m* associated 
with neutral stability was independent of the varying parameters when ft 
was very near zero. Letting ft approach zero in equation (8a) gives the 
value 


m' = 


Y a 0 + V < 


x a 

1 + — 
u 


Setting x a /u equal to zero, as was done for figures 3 and k, reduces the 
expression to 


m* 



The location of the boundary for small n, then, is dependent only on the 
aerodynamic derivatives and the ratio x a /u. For negative values of 
xa/u(> -l) the boundary near 0. = 0 is above m' = Y a /Z^ and for 
positive values it is below, as illustrated :Ln figure°5. ° 
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When x a /u = 0 the portion of the boundary near ft = 0 corresponds 
to the point (Y aQ /Z ao , Y ao /Z aQ ) on the boundary of figure 7- Under this 
condition calculations indicate that the component of pitch in phase with 
bending velocity exactly cancels the angle of attack at the fuselage sta- 
tion (q = 0) resulting from the vertical velocity -m'y/V. Changes in 
x a /u move this point as described earlier. 

Knowledge of the behavior of the boundary near ft = 0 allows qualita- 
tive prediction of the effect of changes in Y aQ /Z ^ in this region, a 
variable not yet considered. 


Minimum Value of m' Necessary for Undamped Oscillation 


All the boundaries of figures 3, 4, and 5, except the one for the 
canard configuration in figure 5, were found to correspond to locations 
on the right-hand boundary in figure 7. It was also mentioned earlier 
that the combination of parameters necessary to reach the boundary for 
the canard configuration in figure 5 was a very likely one, and further 
examination of the stability criterion indicates that almost any combina- 
tion of parameters necessary to achieve neutral stability on the left- 
hand boundary of figure 7 is i* 1 the same category. For airplanes of more 
or less conventional configuration then, the right-hand boundary of fig- 
ure 7 will almost certainly be the critical one. 

If only the right-hand boundary is considered, figure 7 shows that 
undamped oscillations cannot be achieved for values of m' less than Y^, 
since the boundary approaches this value asymptotically as 0y increases 
indefinitely. It should be remembered that this was derived with a 
stability criterion which was useful only for ft <1.0, but the minimums 
of the various boundaries calculated always occurred at values of ft less 
than 1.0, so this minimum value of m 1 is a useful quantity. If m' is 
below this value, instability by the mechanism considered in this study 
is very unlikely. Probable values of Yq for a conventional type airplane 
lie between 0.25 and 0 . 30 . 


CONCLUDING REMARKS 


An oscillatory instability in wing bending coupled with airplane 
pitching was examined by means of a simple two degree of freedom repre- 
sentation of a straight -winged airplane. Although many effects were not 
accounted for, it is felt that the basic mechanism of the instability 
was well represented in the system. The results should be considered a 
guide to more detailed analyses for particular cases because of the sac- 
rifice of detail in this study in the interest of generalization. The 
effects of changes in several quantities describing the airplane and its 
flight condition on wing-bending stability are shown in the form of 
stability boundaries. 
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It was shown that for practical purposes there is a minimum value of 
the ratio of effective wing- tip mass to total airplane mass below which 
wing-bending instability cannot occur by this mechanism. It was also 
found that this instability is most likely to occur under conditions of 
high ratios of effective tip mass to total airplane mass, ratios of 
natural coupled wing-bending frequency to natural uncoupled pitching 
frequency near or below 1.0, and tip mass center-of -gravity locations 
ahead of the total airplane center of gravity. Movement of the tip mass 
center- of -gravity location appeared to be the most effective method of 
avoiding a region of instability for an airplane which is critical in 
this regard. 

Ames Research Center 

National Aeronautics and Space Administration 
Moffett Field, Calif., Sept. 17, 1958 
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Positive 

force 






Figure 1.- Pitch and wing-bending coordinates deflected in a 

positive direction. 
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Coupled bending frequency ratio, £1 


(d) kg = O.lOj u' = 0.?0 
Figure 3 • “ Concluded. 
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(a) Spanwise vertical velocity distribution, ya(r]). 



(b ) Spanwise angle-of -attack distribution, y 


a(ri) . 



(c) Spanwise air load ^ l a 


due to angle-of-attack distribution, 


■ 

y 

v 


a(rj) . 



(d) Spanwise air load due to pitch in phase with bending velocity, 



(e) Effective spanwise air load in phase with y* 

Figure 6.- Contributions of wing-bending velocity and airplane pitch to 
the generalized damping force in the coupled wing-bending mode. 
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Figure 7 *“ Qualitative stability boundary shoving the importance of the 
component of pitch in phase with bending velocity in determination 
of wing-bending stability. 
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